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0O ' ABSTRACT 

On ' 
G\ 

■ Partially phase-mixed structures in galaxies occupy a complex surface of di- 
mension D in six-dimensional phase space. The appearance of such structures to 
observers is determined by their projection into a subspace whose dimensionality 
K is determined by the number of observables (e.g. sky position, distance, radial 
velocity, etc.). We discuss the expected dimensionality of phase-space structures 

\Jq • and suggest that the most prominent features in surveys with K > D will be 

stable singularities (catastrophes). The simplest of these are the shells seen in 

■ the outer parts of elliptical galaxies. 
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1. Introduction 



Of 

The evolution of the phase-space density of stars in galaxies is determined by the collisionless 



Boltzmann equation, which states that phase-space flow is incompressible. Thus a cloud of 
stars in phase space becomes more and more distorted as it evolves; the local or fine-grained 



density around any point in phase space remains the same, but the coarse-grained density 
evolves towards a stationary state. This process, known as phase mixing, is the principal 
mechanism by which stellar systems reach coarse-grained equilibrium. The standard analogy 
is stirring a glass containing 20% rum and 80% Coke; eventually every finite volume element 
in the glass contains 20% rum and 80% Coke, even though infinitesimal volume elements are 
either 100% rum or 100% Coke (e.g. [Arnold and Avez 1968| ). 



Galaxies are not thoroughly phase-mixed, since they are at most a few dynamical times 
old in their outer parts, and since mergers and star formation continuously add new stars. 
There is a variety of direct observational evidence for incomplete phase mixing: (i) Sharp- 
edged features ("shells") in the outer parts of at least 10% of elliptical galaxies are believed 
to arise from the recent tidal disruption of small galaxies (e.g. Hernquist and Quinn 1988| ). 
(ii) Proper-motion and radial-velocity surveys of the local halo reveal clumping in phase 
space flMajewski et al. 1996| ). (iii) The tidally disrupted Sagittarius dwarf galaxy ( |lbata et| 



al. 1997| ) provides a concrete example of a cloud in the early stages of phase mixing, (iv) 
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Moving groups in the solar neighborhood may result from dissolution and mixing of star 
clusters and associations ( [Eggen 1965| , pehnen 1998| ). 

Incomplete phase mixing is likely to play a growing role in the interpretation of observa- 
tions, for several reasons: radial-velocity and proper-motion surveys are rapidly improving in 
size and quality; large telescopes with improved image quality will permit us to examine the 
surface-brightness structure of elliptical galaxies with high spatial resolution and signal-to- 
noise ratio; the Hipparcos mission has dramatically improved the precision of measurements 
of the phase-space distribution of stars in the solar neighborhood, and future space-based 
astrometric missions such as SIM and GAIA will provide proper-motion and parallax mea- 
surements over much of the Galaxy; and the Sloan Digital Sky Survey will soon provide a 
stellar database of unprecedented size and uniformity (10 s stars with five-color photometry). 

Virtually all theoretical descriptions of phase mixing have focused on mixing in a two- 
dimensional phase space, which is the simplest to visualize (Fig. |l|). This paper discusses 
some of the geometrical features of partially phase-mixed systems (we shall call these "phase 
structures") with higher dimensionality, and the projection of this geometry into the observ- 
ables that are measured in a survey. We shall argue that this geometry can be organized by 
two integers: the dimension of the phase structure and the dimension of the survey. 

2. The dimension of phase structures 

A cloud of stars orbiting in an integrable potential (the "host galaxy") can be viewed as a 
surface or manifold of dimension D < N embedded in iV-dimensional phase space (usually 
N = 6). For example, a cold disk has D = 2 (the disk occupies a two-dimensional surface 
in configuration space, and the velocity at each position is unique), collisionless cold dark 
matter has D = 3 (three spatial dimensions but zero random velocity), and a hot galaxy has 
D = 6. Since real stellar systems always have non-zero thickness and velocity dispersion, we 
are referring to an "effective" or "coarse-grained" dimension; loosely speaking, a structure 
in a phase space of iV dimensions has effective dimension D if its extent in D independent 
directions is much larger than its extent in the other N — D directions. 

Because the flow of the cloud through phase space is incompressible, the dimension 
of the manifold does not change in the initial stages of phase mixing, although the shape 
of the manifold becomes more and more complicated. However, the effective dimension 
can change over longer times, through several distinct processes. The most obvious is that 
the phase-mixing scale grows smaller and smaller with time (see Figure |l]), and eventually 
becomes smaller than the coarse-graining scale. At this point the phase structure has the 
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Fig. 1. — Phase mixing in a two-dimensional phase space. The figure shows the evolution of 
5000 points following the equation of motion x = —y/r, y = x/r where r 2 = x 2 + y 2 . 
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same (coarse-grained) dimension as the phase space, unless all of the stars in the cloud have 
the same value of one or more isolating integrals of the motion. Gravitational scattering 
by small-scale irregularities (e.g. massive objects) can also lead to diffusion of the phase 
structure and an increase in effective dimension. 

Less obviously, phase mixing can also decrease the number of effective dimensions. For 
example, consider the cloud resulting from the disruption of a small hot galaxy. The motion 
of the stars in the cloud can be described by a Hamiltonian H(J), where (w, J) are action- 
angle variables in the host potential. The equations of motion are 

dw OH (J) . . 

J = const > *- = -ar sn(J) - (1) 

The trajectory of the centroid of the cloud is (J c (t), w c (i)) = (Jo, w + n t), where O = 
fi(Jo) and (Jo, w ) is approximately the location of the center of the satellite at the disruption 
time t — 0. 

Similarly, the trajectory of a star in the cloud is (Jo + A J, w + Aw + O(J + AJ)£), 
where |AJ| <C J and |Aw| <C 1. This can be simplified to 

(J(t), w(t)) = (J c , w c (t)) + (AJ, Aw + HAJt) + 0(AJ, Aw) 2 , (2) 

where 

H lk = d 2 H(J )/dJ l dJ k (3) 

is the Hessian of the Hamiltonian. At large times the extent of the cloud is dominated by 
the terms oc t in equation (0), so we may write approximately 

J(t) ~ J + O(AJ), w(t) ~ w c (t) + HAJt + 0(Aw). (4) 

Since H is symmetric, it is diagonalizable, that is, there exists an orthogonal matrix A such 
that 

AHA" 1 = D(A), (5) 

where Dij(X) = Aj5y, A = {A.;} are the eigenvalues of H, and A* = A -1 . We now make a 
canonical transformation to new action-angle variables (J', w') using the mixed generating 
function S(J', w) = J' Aw; thus 

w'=^ = Aw, J = ^± = A t J / , w = AV, J' = AJ. (6) 
a J' aw 

Equation now simplifies to 

m ~ 4 + O(AJ'), w'M ~ w' a (t) + A,Aj;t + O(Aw'). (7) 
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This result shows that the shape of the expanding cloud is largely determined by the 
eigenvalues A of the Hessian H, which are invariant under the orthogonal transformation 
A. If two of the three eigenvalues are zero (say A 2 = A 3 = 0) then the cloud expands only 
along the w[ direction in phase space, yielding a one-dimensional tidal streamer that grows 
linearly with time. If one of the eigenvalues is zero (say A3 = 0) then the cloud expands in 
two dimensions in phase space. However, if in addition |AiAJ{| | A 2 A ^ I or vice versa, then 
the expansion will effectively be one-dimensional. Thus the dimension D has been reduced 
from 6 to 1 or 2, because the cloud expands much faster in some angles than in others, and 
not at all in action space. 

We illustrate these remarks with some examples. The Hamiltonian for the triaxial 
harmonic oscillator may be written H(J) = J2k=i w kJki where is the frequency along axis 
k. The Hessian for this Hamiltonian is zero, so the disrupted system does not expand at 
all: the frequencies are independent of the actions and phase mixing does not occur. The 
Hamiltonian for the Kepler potential is H(J) = — \{GM) 2 j J\ , where J\ = (GMa) 1 / 2 and a is 
the semi-major axis. The eigenvalues of the Hessian are (— 3/a 2 , 0, 0), so the expanded cloud 
is one-dimensional, even if its original state was 6-dimensional. Simulations of disruptions 
of small galaxies in a spherical logarithmic potential yield one-dimensional clouds as well 
( [Johnston et al. 1996| ), but the precise conditions determining the effective dimension of the 
cloud resulting from tidal disruption of a small stellar system have not yet been studied. 



3. Catastrophes 

We now survey the phase space by measuring K of the N phase-space coordinates of each 
star (we call this A-dimensional space the "observable space"[]), and ask how to detect the 
D-dimensional phase structure. For example, an image or set of sky positions has K = 2; a 
data cube (positions and radial velocities) has K = 3; a proper-motion survey has K = 4, 
etc. A bolometric laboratory detector of cold dark matter measures only their energy and 
moreover is restricted to a single point in configuration space. 

If D < K the phase structure will appear as a distinct entity in the observable space, 
which can be detected by standard techniques: plotting the data on paper (K = 2), visu- 
alization software (K = 3), cluster- finding algorithms, etc. More specialized and powerful 
techniques can be used if the phase structure has known properties: for example, the dis- 
ruption of distant satellite galaxies in a spherical potential leads to phase structures which 
are great circles on the sky ( [Lynden-Bell and Lynden-Bell 19951 , [Johnston et al. 1996|) . 



Other names include control space or external variables. 
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In this discussion we concentrate the more challenging case when D = K, so that the 
phase structure covers most or all of the observable space. In this case, the most prominent 
features in the survey will be the singularities of the phase structure, that is, the locations 
where the projection of the phase structure into the observable space leads to a singularity. 
For example, if in Figure [l] the observations yield only the horizontal coordinate, then near- 
singularities arise in the plot at t = 50 when the curve representing the phase structure is 
vertical. 

We can generalize this simple concept to singularities with higher dimension usins] catas- 
trophe theory flPoston and Stewart 1978| , [Berry and Upstill 1980| , pilmore 1981| ). We begin 



by setting up a coordinate system . . . ,Ud) on the phase structure that contains the 
cloud. We complete the coordinate system of the phase space by (wd+i, . . . ,un), which are 
chosen so that ud+i = • ■ • = un = on the phase structure. Thus the phase structure is 
specified by 

M^ =0 , j = l,...,N where 0(u) = f) u). (8) 

aU 3 j=D+l 

The advantage of specifying the phase structure in terms of the generating function <fi is that 
is a single-valued function of the phase-space coordinates. 

Next let (xi, . . . ,xd) be the coordinates of the observable space, and complete the 
coordinate system of the phase space by (yi, . . . , Un-d) (we call this (N — D)-dimensional 
space the "hidden space" 0). For given values of the observable coordinates, the location or 
locations of the phase structure are specified implicitly by the equations 

^>=0, ./ ! V II (9) 

Singularities in the observable space arise when there are displacements in the hidden space 
that leave (0) unaltered, that is, for which 

if y? X) ^* = 0. j = l,...,N-D. (10) 



k=l 



dyjdy k 



These equations have a non-trivial solution if the determinant of the (N — D) x (N — D) 
Hessian matrix vanishes, that is, if 

det^>=0. (11) 



2 Other names include state space or internal variables. 
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We can restrict our attention to structurally stable singularities (roughly speaking, ones 
for which small changes in lead to small changes in the locus of the singularity); the reason 
is that unstable singularities represent a subset of measure zero and hence are atypical, at 
least for N < 7 QZeeman E?77|). 

Catastrophes are structurally stable singularities of gradient maps. They are classified 
by their codimension, which is the dimension of the observable space minus the dimension 
of the singularity Thus codimension 1 corresponds to a singular point in a survey with 1 
observable, a singular line in a survey with 2 observables, etc. The corank of the catastrophe 
is the minimum value of N — D, the dimensionality of the hidden space, which is required 
for the singularity to occur. A third classification is by the degree k of the singularity: the 
mean density in a neighborhood of the observable space of radius e around the singularity 
is oc e~ k . 



3.1. Fold 

The fold is the only catastrophe with codimension 1. Since the corank of the fold is also 1, the 
simplest example occurs in a phase space with two dimensions, one of which is observable. 
For example, the generating function 

4>(y; x) = ly 3 + xy (12) 

implies that the phase structure is the parabolic curve x + y 2 = 0, and this manifold has 
a fold catastrophe at x — 0. If the linear density on the curve is uniform, the observable 
density is 

p(x) cx (l - —J , x < 0, (13) 

and zero otherwise. The square- root divergence in density on one side of the singularity is 
characteristic of a fold, and implies that the degree k = |. 

Examples of fold catastrophes occur in the projection of a spherical shell onto a plane, 
rainbows, sunlight sparkling on the sea, twinkling of starlight, gravitational lensing, etc. 

The simplest examples of fold catastrophes in phase mixing are shell structures in el- 
liptical galaxies, which are believed to arise from disrupted companion galaxies. Shells are 
one-dimensional structures that appear in two-dimensional observable space (the two coordi- 
nates on the sky plane) and hence if they are singularities they must have codimension 1 and 
hence must be folds. Hernquist and Quinn (1988) carried out simulations of the evolution 
of disrupted companions in spherical host galaxy potentials, and several of their conclusions 
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can be interpreted in terms of the results we have derived so far. They find that shell forma- 
tion requires the accretion of either a dynamically cold or a spatially compact companion; 
these are precisely the conditions required so that the cloud dimension D < 2 (§0), which 
in turn is required so that singularities are present in a two-dimensional observable space. 
They distinguish "spatial folding" (projection from three-dimensional configuration space 
to two-dimensional space) from "phase-wrapping" (projection from phase space into three- 
dimensional space) but this distinction is not fundamental: both are projections from phase 
space into observable space. Finally, they point out that some shells arise from step functions 
in surface density rather than square-root singularities; the former are not singularities and 
hence are not described by catastrophe theory. 



3.2. Cusp 

The cusp is the only catastrophe with codimension 2, and has corank 1. The cusp is a singular 
point in a 2-dimensional observable space, a line in 3-dimensional observable space, etc. For 
example, consider a phase space with N = 3 dimensions, two of which are observable. The 
generating function 

0(y; x ) = ;b 4 + \ x iv 2 + x iv ( 14 ) 

implies that the phase structure is given by the surface y 3 + x 2 y + x\ = 0, which is singular 
(in the sense of eq. |ll|) along the lines 

Xl = ± ^ (- x 2) 3/2 • (15) 
Each of these lines is a fold catastrophe, and their junction is a cusp catastrophe. 

If the surface density of stars is uniform on the phase structure, the surface density in 
the observable space (x%,x 2 ) i s gi ven parametrically by 

1/2 

, where x\ = —y 3 — x 2 y. (16) 

At a distance Ax 2 from the fold (|Ax 2 | "C \x 2 \ ~ |xi| 2 / 3 <C 1), the density £ oc | | ~ 2 / 3 for 
Ax2 > and £ oc |xi| _1 / 3 (— Ax 2 ) -1 ^ 2 for Ax 2 < (the characteristic fold behavior). The 
degree of the cusp singularity is k = |. 

To illustrate the appearance of a cusp catastrophe, the first three panels of Figure |] 
show a survey that samples the surface density (|1^) with 300, 1000, and 3000 data points 
in the interval X\,x 2 6 [—1, 1]. The fourth panel shows a higher-resolution survey with 3000 
points in x\,x 2 G [—0.1, 0.1]. The general cusp in a 2-dimensional observable space is formed 
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Fig. 2. — The cusp catastrophe defined by the generating function (|14|) , assuming a uniform 
density of stars on the manifold defined by d<p/dy = 0. The first three panels show samples 
of 300, 1000, and 3000 points over the square [—1, l] 2 ; the fourth panel shows a sample of 
3000 points over [—0.1, 0.1] 2 . The cusp point at the origin unfolds into two fold catastrophes. 
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from this special case by smooth distortion (i.e. diffeomorphism) of the observable space, 
and still consists of two fold catastrophes meeting with a common tangent, as in equation 

(0)- 

Examples of cusp catastrophes include the caustic curve seen on the surface of coffee in a 
cup, and the critical point of the van der Waals equation of state in the pressure-temperature 
plane. 

One might speculate that some of the dwarf spheroidal satellite galaxies are cusp catas- 
trophes rather than bound equilibrium stellar systems, but their isopleths are approximately 
elliptical ( [Irwin and Hatzidimitriou 1995| ) and do not resemble the density contours that are 
expected near a cusp. 

3.3. Catastrophes with codimension 3 

There are three catastrophes with codimension 3, the swallowtail, elliptic umbilic, and hy- 
perbolic umbilic. The swallowtail has corank 1, and thus first appears in a phase space with 
N > 4 dimensions; the umbilics have corank 2 and thus require N > 5. All of these singular- 
ities are points in a 3-dimensional observable space, lines in 4-dimensional observable space, 
etc. Just as a cusp occurs at the junction of two fold curves, these occur at the junction of 
cusp curves, which are connected by fold surfaces. 

3.3.1. Swallowtail 

Consider a phase space with N = 4 dimensions, three of which are observable. The generating 
function 

<Pi.y^) = \y b + + \x 2 y 2 + x x y, (17) 

implies that the phase structure is given by y A + x 3 y 2 + x 2 y + x\ — 0, which is singular (in 
the sense of eq. |Tl|) along the surfaces defined in terms of the parameters (y,xs) by 

xx = 3y 4 + x 3 y 2 , x 2 = -2{2y 3 + x 3 y). (18) 

Each of these surfaces is a fold catastrophe; the folds meet at cusp lines defined by 

x 1 = -3y\ x 2 = Sy 3 , x 3 = -6y 2 . (19) 

The fold surfaces intersect one another along the parabola x\ = Ax\. The degree of the 
swallowtail is k — §. 

Figure |3| shows the fold surfaces and cusp lines for the swallowtail catastrophe. 
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3.3.2. Elliptic umbilic 

Consider a phase space with N = 5 dimensions, three of which are observable. The generating 
function 

4>{Vu V2\ x) = y\ - 3y x y\ - x s(Vi + vl) - X 2V2 - sij/i (20) 
implies that the phase structure is given by the surface 

3y{ - 3y\ - 2x 3 y 1 - x x = 0, 6yiy 2 + 2x 3 y 2 + x 2 = 0. (21) 

which is singular along the surfaces defined parametrically by 

Xl = 3y 2 1 -3y 2 2 ±6y 1 (y 2 1 +y 2 2 ) 1 ^ J x 2 = -Qy 1 y 2 ±Qy 2 (y 2 1 +y 2 2 )^ 2 , x 3 = T (y 2 + y 2 2 ) 1 / 2 . (22) 

Each of these surfaces is a fold catastrophe; the folds meet at three parabolic cusp lines 
defined by x\ + x\ — x\ and confined to planes at angles of 120 deg that intersect along the 
x 3 -axis (see Figure £|). The degree is k = 1. 



3.3.3. Hyperbolic umbilic 

The generating function 

0(2/i, V2\ x) = y\ + y\ - x 3 y x y 2 - x 2 y 2 - x x y x (23) 
implies that the phase structure is given parametrically by 

X! = 3y\ + ei(2/i2/ 2 3 ) 1/2 , x 2 = 3y 2 2 + ti{yly 2 ) 1/2 , x 3 = e 2 (ym) 1/2 , (24) 

where q = ±1 and y±, y 2 > 0. Each of these surfaces is a fold catastrophe; the folds intersect 
along the coordinate axes x\ > 0, x 2 = x 3 = and x\ = 0, x 2 > 0, x 3 = (Fig. £|). The cusp 
line is defined parametrically by 

xi = x 2 = 9y 2 , x 3 = 6e 2 y, y > 0. (25) 

The degree is k = 1. 

These exhaust the catastrophes with codimension < 3. Catastrophes with higher codi- 
mension are only relevant when the observable space has dimension D > 4. Catastrophes of 
higher codimension become increasingly complicated, and are less important since we have 
argued in §|| that typical phase structures will often have dimension 3 or less. 




Fig. 4. — The elliptic umbilic (left) and hyperbolic umbilic (right) catastrophes. 
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4. Summary 

It is important to understand the geometry of phase structures before searching for them. 
We have focused on two specific aspects of this geometry: the expected dimension of the 
phase structure in (usually 6-dimensional) phase space, and the classification of singularities 
when a .D- dimensional phase structure is projected into a D-dimensional survey. 

Phase mixing of disrupted small, hot galaxies generally leads to 3-dimensional phase 
structures but the dimensionality can be smaller if the potential of the host galaxy has 
special symmetries (e.g. spherical or Keplerian), or if the eigenvalues of the Hessian (eq. ^) 
are very different in magnitude. Phase mixing of large, cold galaxies leads to 2-dimensional 
phase structures. 

The structurally stable singularities in D-dimensional surveys are folds (D > 1), cusps 
(D > 2), swallowtails, elliptic umbilics and hyperbolic umbilics (D > 3). 

Even though phase mixing is a simple process, there remain many unresolved theoretical 
issues. How fast do small-scale gravitational irregularities or large-scale orbital chaos disrupt 
phase structures? Do phase structures provide a significant source of relaxation in galaxies 
( [Trcmaine and Ostriker 1998j )? What are the properties of phase structures in galaxies that 
we expect from standard models of structure formation? What statistical measures can we 
use to characterize phase structures and their projections? 
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